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I . Introduction 

The relative effectiveness of two motion strategies to 
reduce a target's vulnerability to localization is analyzed 
in this report. A target is a patroling submarine and the 
vulnerability arises because of the ability of a surveillance 
system to determine the submarine's position from time to time. 

The analysis is based on the following model: A target 

is exposed to a surveillance system periodically. The distance 
the target moves during an exposure is essentially zero, but 
the exposure is long enough so that detection of the target 
by the surveillance system can occur. The target is aware of 
its exposures, but the exposures are unavoidable. 

When a target is detected by the surveillance system, the 
target's position is determined, but no other information is 
obtained. The position, which will be called a datum, deter- 
mines a search region in a plane that is bounded by circle of 
radius u^^t that is centered on the datum where u^ is the tar- 
get's maximum patrol speed and t is the time since detection. 

A searcher is cued by the surveillance system; but, due to 
various causes, a delay (time late) occurs before the start 
of a search. 

Although its exposure times are known to a target, infor- 
mation about the location and the activity of searchers is 
not. However, the number of searchers and the searchers' 
sweep rates are small enough so that, without a surveillance 
system cue, the searchers do not represent a threat to the 
target . 
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The first motion strategy is a member of a class of motion 
strategies. For this class, a target course and speed are 
determined by independent stochastic processes. At each expo- 
sure time, the target's course 9 is determined by the uniform 
distribution defined by the density function f0(9) = 1/2-n 

where 0 < 9 < 2 tt and its speed u is determined by the triangular 

2 

distribution defined by the density function fy(u) = 2u/u^ 

where 0 < u < u . 

- - m 

The second motion strategy is also a member of a class of 
motion strategies. For this class, a target course is deter- 
mined by a stochastic process but its speed is deterministic 
and remains constant throughout the motion. At each time 
determined by a Poisson process with parameter a, the target's 
course 6 is determined by the same uniform distribution that 
determines the course in the first class of strategies. 

The members of the first class are identified by the para- 
meter u^ and the exposure period. The members of the second 
class are identified by the parameter u, the target's speed, 
and the parameter a. The two classes are defined in more 
detail in Appendices 1 and 2. 
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II. The First and Second Motion Strategies 

Danskin has shown in Reference 1 that the strategy of 
the first class with the parameter equated to the maximum 
patrol speed u,^ is an optimum strategy under certain conditions 
The conditions are those that exist in the model for t < t 

e 

where t is the time from a detected exposure and t^ is the 
time between exposures. During this interval, the strategy 
is optimum in the following sense: Given a target's motion 

strategy is knov/n to a searcher, then, for t < t^ and prior 
to the start of a search, the searcher's information about the 
target's position is a minimum; since, for the searcher, the 
target's position is uniformly distributed over the region 
bounded by the circle of radius u^^t that is centered on the 
datum. All the remaining members of the first class and all 
the members of the second class are inferior to it in this 
sense. This strategy will be called the first strategy. 

Based on Washburn's results in Reference 2, for a searcher 
that knows a target's strategy, the target's position distri- 
bution for a strategy of the second class is uniform over the 
region bounded by the circle of radius ut centered at the 
datum at t given the target has made exactly tv/o course changes 
by that time. This suggests that the motion strategy of the 
second class should be considered with speed u = Uj^ and at^ = 2 
Here t^ could represent an average time between exposures if 
the exposures were not periodic. This strategy will be called 
the second strategy. With oitg — 2, a most probable number of 
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course changes by is 2; and, after a detected exposure, 
if there has not been a subsequent detection for a searcher 
that knows the target's strategy, a most probable position 
distribution at t^ for the target is the same as the position 
distribution at t^ for a target using the first strategy. 
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III. Some Comparisons of the Two Motion Strategies 

To establish the basis for the first comparison 
of the two motion strategies, consider the following localiza- 
tion search: At a time t after a detected exposure, a sensor 

is placed at the datum. The detection range of the sensor is r^, 
and the time the target can be observed is short enough so that 
during that time the target's movement is negligible. This 
search suggests as a measure of motion strategy effectiveness 
the probability p(t) at t that a target is not on or within 
the circle of radius r^ < u^^t that is centered on the sensor; 
and it will be used for the first comparison. 

In terms of F , . (r;t) the cumulative distribution func- 
K ( t ; 

tion for a target's range from a datum, p(t) = 1 - (r^;t). 

Figures 1 and 2 in Appendix 3 give graphical comparisons of 
F , . (r;t) for the first and second strategies at times equal 
to t^, 2t^, 3t^ and 4t^ given no subsequent detections. As 
can be seen from Figures 1 and 2, for the first comparison and 
these times, the second strategy is superior to the first. 

It can be conjectured that this is the case for all t > 0. 
However, since the measure of effectiveness rewards strategies 
that concentrate probability at u^t, the maximum possible 
target range, the following strategy is superior to both stra- 
tegies for r^ < At the beginning of the motion, choose 

any course but choose the speed equal to Uj^; then hold the 
course and this speed throughout the motion. Clearly, the 
measure of effectiveness is not appropriate for every kind 
of localization search. 
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To establish a basis for a second comparison, consider 
a localization search that is identical to the search that 
was considered for the first comparison except for the place- 
ment of the sensor. For this search, the sensor's placement 
is not restricted to the datum. This search also suggests 
as a measure of motion strategy effectiveness the probability 
p(t) at t that the target is not on or within a circle of 
radius r^ < Uj^t that is centered on the sensor. Note, however 
that for this comparison its value will depend on the sensor 
placement strategy of the searcher. To account for this, 
p (t) will be used for the second comparison where p_(t) is 
the minimum value of p(t) as a function of sensor placement 
strategy. 

Consider a bounding circle of radius u^^t - r^, r^ < 

centered on a datum. For the first strategy, the probability 

,p(t) that a target is not on or within the circle of radius 

r^ is independent of the location of its center if the center 

2 2 

is on or within the bounding circle and p(t) = 1 - r^ /(u^^t) 

for t < t^. If the center is outside the bounding cirole, 

2 2 

p(t) > 1 - r^ / (Uj^t) . Therefore, for the first strategy, 

2 2 

the measure of effectiveness p (t) = 1 - r /(u,.t) for t < t 

'^m c M e 

For the second strategy, the evaluation of in 

terms of r and u,.t is not as tractable. However, it is shown 
c M 

2 2 

in Appendix 5 that p(t) < 1 - r^ /(Uj^jt) for certain values 
of t and r^. Therefore, for the second comparison and these 
values, the first strategy is -superior to the second. 
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Suppose as a placement strategy a sensor is placed randomly 

on or within a circle of radius u,,t + r centered on a datum. 

M c 

Then the probability that the sensor will be at or v/ithin a 

distance r^ from a point on or within a circle of radius u.^t 

2 2 

centered on the datum is r /(u,,t+r ) . In this case, 

c M c 

2 2 

p(t) = 1 - r^ /(Uj^t+r^) regardless of the motion strategy. Con- 
sequently, for any motion strategy, for r^ < u^^t, 

i 1 - r^2/(u^jt+r^)^ For r^ > u^t, p^(t) = 0. 

Both of the comparisons above are related to a search of a 
circular region for a time that is short enough so that during 
the search a target's movement is negligible. A comparison 
related to a search of several circular regions either sequen- 
tially or simultaneously would be more general. In particular, 
a sequential search would allow the consideration of a placement 
strategy in which sensor placement is determined by using the 
information gained as a search progresses. In Reference 1, 
Danskin has analyzed a sequential placement strategy of this 
kind in a game theory context for a target using the first stra- 
tegy in a situation where a datum is established but there are 
no subsequent exposures and consequently no subsequent course 
or speed changes. For t < t^, the analysis provides a means of 
obtaining an upper bound on the probability that a target will be 

contained in one of the circles of radius r at the times and 

c 

places that are determined by a sequential placement strategy. 
However, the analysis does not apply to the second strategy; 
because with that strategy a target's course can change between 
exposures. It appears, as suggested by Danskin in Reference 1, 
that a corresponding analysis for the second strategy would be 
difficult . 
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The model that has been used to analyze the motion 
strategies does not account for factors that could be signi- 
ficant in an operational situation; for example, navigation 
constraints, transit requirements (drift) and evasive maneu- 
vers. In particular, the model does not allow a searcher 
to determine a target's course. 

Suppose the model were extended so that a searcher could 
determine a target's course when the target was at point if 
it detected a target wake at the point. In this case, the 
decay of wake detecability and t^ , the searcher's minimum time 
late, would become critical factors. In particular, suppose 
the decay of at least one kind of wake and t^ were such that 
there was a non-zero probability that a searcher could deter- 
mine a target's course at a datum. Also, suppose the time to 
do this was negligible. Now consider a searcher that had 
determined a target's course at a datum and knew the target's 
motion strategy, but by t < t^ had no other information. With 
the first strategy, the target's position along the target's 
course line would be determined by the triangular distribution 
defined by the density function 

fx(t) = 2x/(Uj^t)^ 

where X(t) is the target's distance from the datum and 
0 < X < With the second strategy, according to Washburn 

in Reference 2, the target's position would be determined by 
the distribution defined by the density function 

^X(t) ,Y (t) = [l/2TT(Uj^t) ^][at/(l-x/Uj^t) ]exp{-at[l- (1-p^) ^] } 
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2 2 3 " 

for p < 1 where p = (x +y ) and the probability exp (-at) 

for the point (ut,0) . The datum is the origin of the coor- 
dinate system and the direction of the positive x-axis equals 
the target's course at the datum. With respect to the first 
and second strategies, for the case being considered, a target's 
choice is between a distribution confined to a line and one 
confined within a circle and at a point. In Reference 3, 

Belkin has argued that of the strategies of the second class, 
the second strategy is also a good choice in this case. 
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IV. Some Conclusions 

The comparisons that are made in Section III emphasize 
the dependence of target motion strategy effectiveness on 
search capability. The game nature of the problem of choos- 
ing an optimum strategy is also emphasized by the second 
comparison . 

To some extent, the choice of the second strategy from 
the second class of motion strategies was arbitrary. Figures 
3, 4, 5 and 6 in Appendix 4 give graphical comparisons of 
^R(t) first strategy and two other members of 

the second class of motion strategies. In Figures 3 and 4, 
the comparison is between the first strategy and the strategy 
of the second class with a = 1/t^. In Figures 5 and 6, the 
comparison is between the first strategy and the strategy of 
the second class with a = 4/t^. As can be seen from Figure 5, 
for the first comparison, the second strategy is superior to 
the strategy with a = 4/t^ for t = t^. For the second com- 
parison, the second strategy is superior to the strategy with 
a = 1/tg based on the arguements in Appendix 5, since for 
this strategy the probability that a target is on the circle 
of radius Uj^t that is centered on a datum is exp(-t/t^) rather 
than exp(-2t/t^) which is the case for the second strategy. 

The difference can be seen by comparing the plots for t = t^ 
in Figure 1 and Figure 3. 

It is not appropriate here to discuss the relevance of 
the model on which the comparisons are based; however, to the 
extent it is relevant, the comparisons provide some basis for 
the choice of a motion strategy. 
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Appendix 1 



"THE FIRST CLASS OF MOTION STRATEGIES" 

The first class of motion strategies is defined as follov/s: 

At the beginning of the motion, choose a course 9^^ from a uniform 

distribution with density function fQ(9) = 1/2 tt where 0 < 9 < 2tt 

and choose a speed u^ from a triangular distribution with density 

2 

function f„(u) = 2u/u where 0 < u < u . Maintain the course 0, 

U m - - m 1 

and speed u^^ until the first exposure. After the first exposure, 

choose a new course Q2 ^ speed U 2 by repeating the above 

procedure. Continue in this manner until the end of the motion. 

Suppose a surveillance system detects an exposure of a target 

using a motion strategy from this class and suppose the strategy 

is known to a searcher. In this case, the information about the 

target's position that is available to the searcher between the 

time of detection and the next exposure is that the target's 

position is uniformly distributed over a plane region bounded by 

a circle with radius u t that is centered on the datum. 

m 

The above statement can be established as follows; Let the 

datum be the origin of a rectangular coordinate system with east- 

west coordinate x and north-south coordinate y. Then a point's 

bearing d and range r from the datum are related to the point's 

rectangular coordinates by x = r sin d and y = r cos d* If the 

point represents the target's position, then relative to the 

searcher all of these quantities are random variables. And, 

for 0 < t < t^, X(t) = R(t) sin d and Y(t) = R(t) cos d where 

R(t) = ut and d = 9. Since U and 0 are independent for 0 < t < t^ 

where t is the time of the next exoosure, their joint density 
e 
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2 

function is given bv f„ p.(u,0) = u/tt(u ) where 0 < u < u 

■ u , y m - - m 

and 0 < 9 < 2tt . Therefore, for the 0 < t < t^ , the joint density 

2 

function of R(t) and 'J' is given by f_,.> ^(r,(f);t) = r/TT(u t) 

K ( tl ; ^ ^ n 

where 0 < r < u t and 0 < 4> < 2tt. And, for 0 < t < t , the joint 
--m - -*e 

density function of X(t) and Y(t) is given by f^Ct) Y(t) ” 

1/7t(u t) ^ where + y^ < (u t)^. 

For 0 < t < t^, the marginal density function for the taraet's 

2 

range from the datum is f„,. , (r;t) = 2r/ (u t) where r < u t. The 
^ R(t) ' m - m 

cumulative distribution function for R(t) is given by 

^R(t) ''''■*=> ' 

where r < u t. 

m 

After a detected exposure, if a target is not subsequently 
detected, the target's position distribution is no longer uniform 
after its next exposure. Histograms which correspond to cum.ula- 
tive distribution functions for the target's range from the ori- 
gin at subsequent exposure times are shown in Appendix 3 . The 
distribution of the target's bearing from the origin remains 
uniform. In the limit, the target's position distribution ap- 
proaches a circular normal distribution centered on the origin. 
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Appendix 2 



"THE SECOND CL/iSS OF MOTION STRATEGIES" 

The second class of motion strategies is defined as follows: 
At the beginning of the motion, choose a patrol speed u. This 
speed will be maintained throughout the motion. Next, choose a 
course 0^ from a uniform distribution with density function 
fQ(0) = 1/2tt Inhere 0 < 0 < 2tt and choose a time t^ from an expo- 
nential distribution with density function f^(t) = a exp (-at) 
for 0 < t. Maintain the course 0^ for a time t^. At the end of 
that time, determine a new course 02 and time t 2 by repeating the 
above procedure. Continue in this manner until the end of the 
motion. 

Suppose a surveillance system detects an exposure of a target 
using a motion strategy from this class and suppose the strategy 
is known to a searcher. In this case, the information about the 
target's position that is available to the searcher at a time 
t > 0 is that the target's position is described by a joint den- 
sity function determined by Washburn in Reference 2. 

The marginal density function for the target's range from the 
datum is ^R(t) ~ ( 1/ut) [atp/ ( 1-p ) '1exp{-at[l- (1-p ) ^1} ; 

where r < ut and p = r/ut. For r < ut, the cumulative distribu- 
tion function for R(t) is given by 

FR(t) (r't) = 1 “ exp{-at[l- (1-p^) } 

The probability the target's range equals ut is exp (-at) . For 
all r, the target's bearing is uniformly distributed between 0 
and 2tt. 
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Plots of 



the target's 
distribution 



^R(t) shown in Appendix 
position distribution approaches 
centered on the origin. 



3. In the limit, 
a circular normal 
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Appendix 3 

"A GRAPHICAL COMPARISON OF THE FIRST & SECOND MOTION STRATEGIES" 

Range cumulative distribution functions for a target's 

range from a datum given no subsequent detections for the first 

and second strategies are shown in Figures 1 and 2. The range 

units are u,.t and the distributions are for times equal to 
Me 

t , 2t , 3t and 4t . 
e' e e e 

The cumulative distribution functions for the first stra- 
tegy are represented by histograms. The histograms v/ere gen- 
erated by using a computer simulation. The cumulative distri- 
bution functions for the second strategy are represented by 
plots. The plots were generated using a computer plot routine 
and the cumulative distribution function given in Appendix 2 
with the course change rate a = 2/t^. 

In Appendix 6, an analytical comparison of the tx^?o range 
distributions is given in terms of the first and second moments 
about the origin. 
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Figure 1. Range cumulative distribution functions at 
the indicated times. The histograms are 
for the first strategy and the plots are 
for the strategy of the second class with 
a = 2/t^ (the second strategy) . 
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Figure 2. Range cumulative distribution functions at 
the indicated times. The histograms are 
for the first strategy and the plots are 
for the strategy of the second class with 
a = 2/t^ (the second strategy) . 
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Appendix 4 

"A GRAPHICAL COMPARISON OF THE FIRST & OTHER MOTION STRATEGIES" 

Range cumulative distribution functions for a target's 
range from a datum given no subsequent detections for the first 
strategy and two members of the second class of motion stra- 
tegies are shown in Figures 3, 4, 5 and 6. 

Figures 3 and 4 show the cumulative distribution functions 
for the first strategy and the strategy of the second class 
with a = 1/t^. Figures 5 and 6 show the cumulative distribu- 
tion functions for the first strategy and the strategy of the 

second class with a = 4/t . 

e 

The range units are and the distributions are for 

times equal to t^ , 2t^, 3t^ and 4t^. The histograms and plots 
were generated in the same manner as those in Appendix 3. 
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Figure 3. Range cumulative distribution functions at 
the indicated times. The histograms are 
for the first strategy and the plots are 
for the strategv of the second class with 
a = 1/t^. 
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Figure 4. Range cumulative distribution functions at 

the indicated times. The histograms are 

for the first strategy and the plots are 

for the strategy of the second class with 

a = 1/t . 
e 
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I 2 

Figure 5. Range cumulative distribution functions at 
the indicated times. The histograms are 
for the first strategy and the plots are 
for the strategy of the second class with 
a = 4/t^. 
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Figure 6, Range cumulative distribution functions at 

the indicated times. The histograms are 

for the first strategy and the plots are 

for the strategy of the second class with 

a = 4/t . 
e 
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Appendix 5 

"A BOUND ON THE SECOND COMPARISON MEASURE OF EFFECTIVENESS" 

A bound on the second comparison measure of effectiveness 

for the second motion strategy can be computed as follows: 

The probability that a target is on an arc length s of a circle 

of radius Uj^t that is centered on a datum is {s/ 27 rUj^t) exp ( -2t/t^) . 

The maximum arc length that can be intercepted by a circle of 

radius r < u..t is equal to s., = 2u.,t sin ~(r /u„t) . Therefore, 
c M ^ MM c M 

the probability that a target will be on or within a circle 
of radius r that is positioned so that it intercepts an arc 
length s^^ is at least equal to 

(1/tt) sin~^ (r^/Uj^t) exp (-2t/tg) . 

2 

For the first strategy, the probability is equal to (r^/Uj^t) . 

Setting x = r^/Uj^t and equating these two expressions, one 
has the following equation in x: 

x^ = (1/tt) (sin ^x) exp (-2t/tg) . 

Since for the second strategy only the probability that the 
target is on an arc is being considered, the solutions to this 
equation are such that at least for r^ < xu^^t, the first motion 
strategy is superior to the second. For t = t^ , x = .0431; 
for t = t /2, X = .1174 and for t = t /4, x = .1943. 

G ^ 
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Appendix 6 

"AN ANALYTICAL COMPARISON OF THE FIRST & SECOND MOTION STRATEGIES" 

An analytical comparison of the range cumulative distribution 
functions for a target's range from a datum given no subsequent 
detections for the first and second strategies is given in this 
appendix. The comparison is in terms of the first moment mj^(t) 
and second moment m 2 (t). Because of the nature of the distribu- 
tion for the second strategy, it is convenient to use the follow- 
ing formulas: 

mj_(t) = / [l - (r;t) ]dr 

0 






m2 ( t) 



= / 2r[l - 






For the strategies of the first class, for t < t^, 

m, (t) = 2(u t)/3 and m- (t) = (u t)^/2. 

1 m z m 

For the strategies of the second class. 



m 



tt/2 

, (t) = (ut) exp(-at) / exp (at cos9) cos 6 d0 
^ 0 



This can be shown by first substituting p = r/ut, then sin9 = p. 
The integral could not be evaluated in closed form; however, it 
was evaluated numerically for specified values of a and t. In 
addition. 



m2(t) = 



(u^/a^) [exp (-at) - 1 + at] 



This can be shown by first substituting p = r/ut and then 



2 1 ^2 

X = 1 - p . 
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Evaluating m 
tegies at t = 
strategy, 
second strategy, 
Clearly, tor t = 
distributions do 



(t) and ni 2 (t) for the first and second stra- 

gives the following results. For the first 

2 

= .67{u,,t ) and m_(t ) = .5(u,.t ) . For the 

Me 2 e Me 

t , the first and second moments of the two 
e 

not differ significantly. 
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